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Abstract: We consider some of the recent proposals in which weak-scale dark matter is 
accompanied by a GeV scale dark sector that could produce spectacular lepton-rich events 
at the LHC. Since much of the collider phenomenology is only weakly model dependent it 
is possible to arrive at generic predictions for the discovery potential of future experimental 
searches. We concentrate on the production of dark states through bosons and electroweak- 
inos at the Tevatron or LHC, which are the cleanest channels for probing the dark sector. We 
properly take into account the effects of dark radiation and dark cascades on the formation 
of lepton jets. Finally, we present a concrete definition of a lepton jet and suggest several 
approaches for inclusive experimental searches. 



1. Introduction 

In light of recent astrophysical observations, the authors of [1] have proposed a broad class of 
theories in which the annihilation of ~TeV scale dark matter in the galactic halo accounts for 
the anomalous excess of cosmic ray leptons. While this dark matter is probably inaccessible 
to colliders, it is accompanied by a ~GeV scale dark sector which couples, albeit very weakly, 
to the standard model (SM). The dark sector states are relatively light and can be produced 
at high energy colliders, only to cascade decay through the dark sector and ultimately return 
to the visible sector as electrons, muons, and possibly pions. The outgoing leptons emerge in 
the detector as lepton jets [2], which arc highly collimatcd multiple muons or electrons that 
result from the decay of these highly boosted dark sector states. As observed in ref. [2] and 
later elaborated on in ref. [3], the addition of supersymmetry also bears with it a number 
of interesting and novel collider signatures. Some of the phenomenology is also relevant for 
scenarios with hidden valleys [4] or where DM decays into a GcV scale sector [5], and as 
such may describe a large class of models. Recently, experimental effort in searching for such 
objects has been reported in [6]. 

The aim of the present work is to provide a quantitative study of the collider phenomenol- 
ogy of a typical event involving the GeV scale sector and resulting in lepton jets. We consider 
the minimal scenario whereby the SM sector and the dark sector are coupled only via kinetic 
mixing terms as detailed below. Effectively, dark sector fermions and scalars have a small 
coupling to the boson and the minimal supersymmetric standard model (MSSM) bino, 
while the dark sector gauge boson couples weakly to the SM electromagnetic current. As 
a result of these couplings, a typical event factorizes into three modular stages, which are 
illustrated in Fig. 1: 

• Electroweak Production 

Once bosons and electroweak- inos are produced, they decay into the dark sector via 
the aforementioned coupling. Rare Z^ decays and electroweak-ino pair production are 
the cleanest channels in which to observe the production of dark sector states due to 
the limited hadronic activity. Decays into the dark sector depend weakly on the dark 
sector details. 

• Dark Sector Evolution 

Once the electroweak states have decayed into the dark sector, the resulting dark states 
are highly boosted and cascade decay down to the bottom of the dark sector spectrum. 
A universal feature of these dark matter scenarios is that there is log enhanced soft 

emission of dark gauge bosons, which increases the multiplicity of light dark sector 
states and ultimately yields a greater number of final state leptons. We simulate these 
dark sector showers and discuss their characteristics. 



1 



• Outgoing Lepton Jets 

Cascade decays and soft emissions through the dark sector result in radiated dark gauge 
bosons which return to the visible sector as collimated lepton jets. After studying the 
shape and distribution of simulated lepton jets and taking into account the possible 
dilution and contamination from the decay of the dark bosons into pions, we suggest a 
concrete definition for a lepton jet which can be used in inclusive experimental searches 
for these objects. While much of the phenomenology we consider is independent of 
the details and spectra of the particular dark sector model, this is not the case for 
the bottom of the dark sector spectrum, which may be probed by studying lepton jet 
shapes. 




Figure 1: A schematic illustration of the type of events we consider in this work. The time evolution 
can be divided into three stages: electroweak boson or -ino production and subsequent decay into the 
dark-sector, evolution through the dark sector, and finally the formation of lepton jets, as delineated 
by the dashed boxes. Such events may also include missing energy. 

In section 2 we review how the dark sector couples to the visible sector and discuss the 
production of dark sector states in rare decays at LEP, the Tevatron and the LHC. We also 
consider electroweak- ino pair production at the Tevatron and LHC. In section 3 we consider 
the evolution in the dark sector which includes dark showering and cascade decays in the dark 
sector itself. Section 4 begins with an analysis of the final state leptons and the formation of 
lepton jets and ends with some proposals for experimental searches. Section 5 contains our 
conclusions. 

2. Electroweak Production 

Let us first review how the visible sector and dark sector are coupled. For a detailed treatment, 
see [3]. As in [1], we assume the existence of a new dark gauge group which contains a C/(l) 
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factor that is spontaneously broken at ~GeV. The associated dark gauge boson kineticaUy 
mixes with the SM photon and the vector-boson. We parameterize these couphngs as 

■^gauge mix = '^]biiv-^ 2^2bfj,iiZ^ (2-1) 

= -y,b,.B'^''-le',b^.Wr (2.2) 

where b^i, denotes the field strength for the dark gauge boson and ei^2 and e[ 2 are related 
by the Weinberg angle. In particular, when only e'l is present, we have ei = e'^ cos 6w and 
62 = e'l sin ^i^. The first parameterization (mass basis) is more useful when discussing SM 
processes while the second (gauge basis) comes in handy when considering the MSSM^ . With 
the addition of supersymmetry, the above coupling implies a kinetic mixing between the dark 
bino, b and the MSSM gauginos, B and Ws: 

/:gaugmo mix = -2ie[b^ a" df,B - lie'^b^ a^O^W^ + h.c. (2.3) 

The gauge and gaugino kinetic mixings can both be eliminated by a set of field redefinitions 
which induce the portal to the dark sector which will be relevant to this collider study: 

>Cportal = eib^J^M + ^2Z^Jb + ^l^Jl + 4^3^g (2-4) 

4 = 9dJ2li {iihldi^hi - hidi'hl) + hlai^K^ (2.5) 

i 

Jl = -iV2gaY,<lih\hi (2.6) 

i 

where Jem is the SM electromagnetic current and Jj, and J-^ are the bosonic and fermionic 
components of the dark gauge current. In particular, J;, contains dark scalar and dark fermion 
bilinears, while contains mixed dark scalar-fermion bilinears. 

Since the dark sector scalars and fermions couple to the boson and the MSSM bino, 
they will be produced in rare Z^ decays and electroweak-ino decays, which we now consider 
in detail. 

2.1 Rare Z^ Decays 

The strongest bound on kinetic mixing with a new ~GeV scale vector-boson comes from the 
muon g — 2 ratio [7] , which constrains the photon mixing to be 

< 3 X 10-5 (2.7) 
^ ~ 100 MeV ^ ' 

On the other hand, the branching ratio of Z^ into a pair of dark fermions charged under the 
dark gauge group is fixed by the mixing parameter, 62, which is not directly limited by this 



'^The operator involving W3 is not gauge-invariant with respect to SUhi^). It should be thought of as 
the result of a higher dimensional operator bij,i,tr{H^W'^'^H)/A^, where A is some high-scale. After the higgs 
condenses one obtains an effective mixing between b^u and Ws^^. 
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bound and can be somewhat larger. The branching fraction of to dark fermions is 



BR {Z ^ ff) 



49i Mzo 



(2. 



while the branching fraction into a pair of dark scalars is a quarter of this. This has implica- 
tions for LEP, as well as the Tevatron and LHC: 

LEP: With a nominal value of = Q!em(-^z) and = 10~^ the 17M on-shell Z^ bosons 
produced at LEP [8] will yield about 150 events of dark fermion pairs. While the LEP data 
can at the very least place stronger bounds on the couplings, it is not clear that existing 
searches covered the type of event topologies involved. Collimated leptons coming off of dark 
sector states may have been missed in exclusive searches requiring isolated leptons. 

Tevatron and LHC: The cross-section for production of a pair of dark sector states off 
the Z^ at Tevatron and LHC is shown in Fig. 2. To a good approximation, the interfering 
diagram involving an off-shell dark vector-boson can be neglected since it produces final 
states which are too soft. With several fb^^ of data, the Tevatron can already probe this 
production channel for BK{Z^ — )• dark sector) > 10^^. The LHC should be sensitive to the 
same parameter region with several hundreds pb^^ of data. 
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Figure 2: The cross-section for production of dark sector states via Z° as a function of the branching 
ratio for both the Tevatron and LHC. The vertical green dashed lines mark the branching ratio for 
£2 = 10"^, 10~2, and 10"^ from left to right, using Eq.(2.8) with ad = ^em = 1/127. 



2.2 Electroweak-ino Pair Production 

Electroweak-ino pair production is well understood and has been discussed extensively in the 
literature. In particular, the tri-lepton final state with its low SM background offers one of 
the better discovery channels for the conventional MSSM [9]. 

In order to minimize the hadronic activity it is useful to focus on the production of 
colorless particles. While strong production channels yield much larger cross-sections, they 
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result in more complicated events. At this early stage of the investigation, we prefer to 
concentrate on the lepton jets by themselves. At the Tevatron and LHC, the colliding particles 
certainly carry color, so the possibilities are few. One possibility is the Drell-Yan production 
of sleptons. Another is the pair production of electroweak-inos, which are superpositions of 
winos, binos, and higgsinos. After these -inos are produced, they will promptly cascade decay 
through the visible sector until they reach the lightest superpartner of the MSSM, which in 
turn decays into dark sector superpartners. The precise details of the MSSM cascade depend 
on the exact MSSM spectrum^, but we will assume that it only produces isolated leptons, as 
is the case for example in decays of the chargino to the neutralino. 

We divide the electroweak production channels into three categories, neutralino-pair, 
chargino-pair, and neutralino-chargino associate production. In Figs. 3 and 4 we depict 
the production cross-sections for the different pure gaugino states at the Tevaton and LHC, 
respectively. The purpose of these graphs is to give an estimate for the cross-sections involved 
for general MSSM spectra rather than concentrate on a particular benchmark scenario. In the 
case of neutralino pair production, the pure bino and wino states are produced only through 
a t-channel exchange of a squark, and their production cross-section is therefore suppressed 
compared with the higgsino state which enjoys its coupling to the Z^. Charged wino pair 
production has a much larger cross-section due to the large coupling of charged winos to the 
electroweak vector-bosons. For the purpose of computing the neutralino-chargino associate 
production we assume a degeneracy between the charged and neutral state, however, in 
realistic spectra those are usually split and the cross-sections are somewhat modified. 

Both Tevatron and early LHC data should be sensitive to large parts of the parameter 
space. It is clear that generically, -ino production is dominated by channels with at least one 
chargino. This fact has important implications on the type of event topologies we can expect 
to encounter in association with lepton jets. As the chargino cascades down to the lightest 
neutralino it will emit a hard isolated lepton or other SM particles as we now discuss. 

2.3 Neutralino Decays into the Dark Sector 

Once produced, electroweak-inos will promptly cascade down to the lightest neutralino, A^i. 
In the process, they may emit leptons, quarks, or Z^/W^ (on-shell or off-shell) depending 
on the precise MSSM mass spectrum. Therefore, some of the events contain isolated leptons 
in addition to the lepton jets generated at a later stage of the event. Such isolated leptons 
should not cause any complications in a properly inclusive search. The cascades to A^i may 
also result in colored particles and hence QCD-jets, but for the purpose of the current study 
we will assume that a substantial fraction of such cascades result in no colored particles. This 
is not a strong assumption as it is satisfied in many concrete examples of MSSM spectra [11], 
and may even be discarded altogether in actual lepton jet searches including QCD-jets. 

^In this work, we will limit ourselves to no more than a single step by assuming that the spectrum contains 
at most one additional ino below the ones produced. It is important to realize that since DM is no longer 
associated with the MSSM per say, any of the MSSM sparticles can be the lightest as long as it is unstable 
and can decay into the dark sector. 
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Figure 3: Production cross-sections for the different ino states at the Tevatron. The left pane includes 
neutralino pair production for the different gauginos. The right pane shows chargino pair production 
as well as neutralino-chargino associated production. A squark mass of 750 GeV was assumed. Cross- 
sections were computed with Pythia [10]. 
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Figure 4: Same as Fig. 3, but for the LHC with center of mass energy of 14 TeV as well as 10 TeV. 

The couplings in eq. 2.4 imply that the lightest neutralino decays into the dark sector via 
its bino or wino fraction. The higgsino state is rendered unstable through its mixing with the 
wino and bino, as depicted in Fig. 5. The complex scalar fields shown in this figure are dark 
gauge eigenstates separated into real and imaginary parts. The precise linear combinations 
of fields which define the mass eigenstates are described in the appendix. In particular, one 
combination is the eaten goldstone boson of the dark gauge symmetry. As we will see, it 
is conceptually simpler to work in the gauge basis, especially when describing the ensuing 
radiation. 

The event shape is determined to a large degree by production, since the electroweak-inos 
are produced almost at threshold. Therefore, the lightest neutralino A'^i is not very boosted 
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Figure 5: The MSSM neutralino can decay to the hght elements in the dark sector. We distinguish 
between the real and imaginary part of the dark higgses, because one linear combination is in fact the 
goldstone boson eaten by the broken dark gauge group. Thus, it can decay directly into lepton pairs 
while the other higgses cannot. 

in the lab frame and so its decay products are not collimated. Each of the two neutralinos 
gives off a pair which yields a total of four highly-boosted dark sector particles which are 
well-separated in the lab frame. They form the seeds of the dark showers described in the 
next section. The final outgoing lepton jets will consist of some isolated leptons coming from 
the electroweak-ino cascades together with up to four lepton jets created by the boosted dark 
sector states. 

We conclude this section by considering some special cases in which the decay into the 
dark sector proceeds slightly differently, resulting in alternative final configurations. For 
instance, in low-scale SUSY breaking scenarios, the gravitino is light and so the lightest neu- 
tralino can decay into a gravitino with the emission of a real photon. This process competes 
with neutralino decays into the dark sector. Generally, those two decay rates are vastly dif- 
ferent [3], but in some instances they may be comparable [12]. If so, the events will contain 
a single hard photon, lepton jets and possibly other isolated leptons. This topology was the 
basis for a recent analysis looking for lepton jets in the Tevatron [13]. 

Second, there are cases where the degeneracy between the charged and neutral winos is 
lifted only by electromagnetic loop corrections, as is the case for instance in anomaly medi- 
ated supersymmetry breaking. The decay — )• W'^ + vr^ proceeds with an approximate 
lifetime of 10"^*^ sec and results in a displaced vertex. Thus, the lepton jets that emerge 
out of the decay of into the dark sector reconstruct this displaced vertex and are them- 
selves displaced. We do not consider such scenarios further in this work, but only comment 
that lepton jets with displaced vertices may require special consideration when designing an 
experimental search because triggers usually assume the leptons to come from the primary 
interaction point^. 

3. Dark Sector Dynamics 

In this section we discuss the evolution of dark sector states that arise from visible sector 
decays. As observed in [3] these states can seed long cascades if the dark sector happens to 

^We thank Roger Moore for pointing this out to us. 
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have a rich particle spectrum. Such a scenario is generic if the dark sector has non-abehan 
gauge dynamics. Along and at the end of these cascades, numerous leptons are produced from 
the resulting dark gauge bosons, yielding highly collimated lepton jets. A primary observation 
of this paper is that dark sectors with relatively sparse spectra, for example abelian theories, 
also provide a sizeable production of lepton jets due to soft particle emission. 

To sec this, consider the scenario described in the previous section, where each neutralino 
arising from a supersymmetric cascade in the visible sector decays into a well-separated and 
boosted dark scalar and fermion pair. Each of these carries an energy of roughly half the 
neutralino mass (~ lOO's GeV). Furthermore, since these dark sector states are so highly 
boosted (and have mass < GeV) they will radiate soft and collinear dark photons with 
high probability [14, 15]. Each dark photon brehmsstrahlung eventually decays into a pair of 
leptons, and contributes substantially to the final lepton multiplicity of the resulting lepton 
jet. This part of the evolution is universal and depends only very weakly on the precise details 
of the symmetry breaking. 

Once the shower completes, we dress the dark scalar s and fermions with the proper mass 
eigenstates as detailed below in section 3.2. Depending on the final state, there may or may 
not be a pair of leptons at the termination point of the dark sector cascade. For example, if 
the final state is the eaten goldstone boson of the dark symmetry breaking, then it will decay 
to a lepton pair directly. If instead it is a scalar which lacks a direct coupling to leptons, then 
it decays through an intermediate (on-shell or off-shell) dark vector-boson (see Fig. 8). 

In this section we discuss each stage of the dark sector evolution: showering, dressing, and 
their combined effect on the structure of the lepton jet. For our analysis we employ the simple 
abelian model described in the appendix with the necessary field content to spontaneously 
break the dark gauge symmetry. As such, much of the abelian model phenomenology is 
indicative of a broad class of models, including more complicated extensions such as non- 
abelian theories. 

3.1 Dark Showering 

Decays into the dark sector produce highly boosted dark particles which radiate soft dark 
photons. The number of photons radiated off a dark higgs, h, or dark higgsino, h, can be 
understood parametrically in terms of the Sudakov double logarithm: 



which is the expected number of soft dark photon emissions occurring within an energy 
window defined by Mdark and Mew- More precisely. Mew is the invariant mass of the 
initial dark sector state, for example the mass of the Z*^ or Ni that decays into the dark 
sector, and M^ark is the the scale of dark sector masses, which regulates the soft and collinear 
divergences. For the estimate in the final step of Eq. (3.1), we have taken Msw/Mdark ~ 10^- 
For M(jark ~ GeV, the only parameters that determine the amount of showering are therefore 
ad and M^^ . Since the ratio of scales. Mew ^ M^ark; is independent of the details of the 




(3.1) 
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GeV-scale spectrum, showering is a universal process. Furthermore, due to the logarithmic 
enhancement, it is clear that showering is an 0(1) effect, even for weakly coupled ^7(1)^. 

In order to perform a detailed study of dark sector evolution, we have implemented the 
virtuality ordered parton shower employed by Sherpa [16] and Pythia [10]. This algorithm 
is given by the repeated application of the following: 1) given an on-shell "mother" particle, 
sample the Sudakov form factor to obtain the mother's new virtuality t and the energy fraction 
z (and 1 — z) going to each of its daughters, and 2) re-shuffle kinematics to be consistent with 
four- momentum conservation. We use the kinematic expressions of [17]. For a substantially 
more detailed description of the virtuality ordered parton shower, refer to the above references. 

For a weakly coupled dark sector with a GeV-scale mass gap, we have found several 
simplifications to the parton shower which have little effect on lepton jet observables. First, 
we note that the dark shower is dominated by splitting functions with both a soft and coUinear 
divergence, h — > hj' and h — > hj'. Other splittings are only enhanced by the coUinear log, 
and therefore constitute a < 10% effect. Second, massive splitting functions and the precise 
GeV-scale virtuality cutoff have little effect. This is because the amount of radiation depends 
only logarithmically on M^ark- We have found that varying the virtuality cutoff of the shower 
by an order of magnitude around a GeV only has a noticeable effect on the pt distribution 
of radiated dark photons for pr GeV, where the dependence is 0(1). Third, the running 
of the dark gauge coupling is negligible, since the theory is weakly coupled in the relatively 
small range between Mew and Mdark) and also constitutes a < 10% effect on the amount of 
showering. Based on the above discussion, we adopt several simplifications for what follows, 
allowing us to keep the showers as model independent as possible. In particular, we only 
include splitting functions that are double log enhanced, we neglect massive splitting functions 
and fix the virtuality cutoff to 1 GeV, and we neglect the running of the dark gauge coupling. 

Lastly, we mention a notable difference between abelian and non-abelian theories per- 
taining to angular ordering. In particular, consider a soft emission of the form A — >■ BC. 
In the case that B and C are both charged (for instance in a non-abelian theory), then any 
subsequent soft emission at large transverse wavelength from either B or C is suppressed by 
a Chudakov-like effect. This is due to interference in the matrix element, and is simulated in 
virtuality ordered parton showers by enforcing 'angular ordering,' where subsequent emissions 
are required to occur at smaller openings angles than prior emissions. In our case, however, 
the theory is abelian, and the dominant process consists of a hard dark-charged h or h line 
emitting dark photons. Since dark photons are neutral, there is no suppression arising from 
angular ordering. 

In Fig. 6 we present the average number of radiated dark photons in rare decays and 
neutralino pair production events. The amount of radiation increases with the dark coupling 
as can be expected. It also increases with larger neutralino mass since the initial dark higgs 
and higgsinos are more energetic. In Fig. 7 we plot the px distributions for the radiated 
photons for different parameters. 
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Figure 6: On the left we depict the normalized distributions for the expected number of radiated 
dark photons in rare decays into the dark sector. On the right is a contour plot of the number 
of soft dark photon emissions per neutralino as a function of the dark gauge coupling ad and the 
neutralino mass A/^. The linear dependence on and logarithmic dependence on Af^ is in accord 
with the naive expectation, Eq. (3.1). The plots were produced using a 3 GeV px cut on the dark 
photons. Both plots are for LHC at 10 TeV center of mass energy. 
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Figure 7: The px distribution of radiated dark photons for rare Z'^ decays and neutralino cascades. 



Altogether, the effects of radiation in rare decays into the dark sector are rather mild. 
The energy scale involved is somewhat lower as compared with neutralino pair production and 
subsequently the radiation is softer and less pronounced. On the other hand, for neutralino 
cascades into the dark sector the effects of radiation are important. The presence of these 
extra radiated photons modifies and enriches the structure of the resulting lepton jets as we 
discuss in the next section. 
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3.2 Decay of dark states back into leptons and pions 

Once the virtuality reaches the dark state's mass (~ GeV) no further radiation is possible and 
the state is placed on-shell. At this point, it may decay back into light standard model particles 
such as leptons and pions^ via the kinetic mixing operators, Eq. (2.1). The characteristics 
of the decay depend on the identity of the dark particle involved since not all are directly 
mixed with the standard model. For illustrative purposes we consider the two dark scalar 
model presented in the appendix, in which case the bosonic spectrum consists of a heavy dark 
higgs, Hd, a dark pseudo-scalar a^, a dark vector-boson and a lighter dark higgs hd- The 
fermionic spectrum might involve the superpartners of these particles, the dark higgsinos and 
gauginos. 

Since only 6^ mixes with hypercharge, it alone can decay directly into leptons. Ha on the 

other hand we take to first decay into two on-shell 6^'s which later also decay, resulting in 4 
leptons. hd can also decay through two 6^'s, but if it is lighter than the dark vector-boson 
then the decay is off-shell and results in a very long lifetime [3, 18]. This typically means that 
hd escapes detection and is counted as missing energy. The pseudo-scalar, ad, can typically 
decay into hd and two leptons through either on-shell or off-shell b^. These different decay 
modes are shown in Fig. 8. While the specifics of the decay are fairly model dependent, we 
believe that the resulting phenomenology it describes is fairly universal in that each dark 
particle can do one of three things: 1) decay into an even multiplicity of leptons (2,4,6 . . . 
leptons); 2) constitute missing energy; 3) a combination of the two. 

In the next section we analyze the structure of Icpton jets that results from the decay of 
the dark states back into leptons. The lepton jets are produced by first simulating the hard 
production process, followed by the decay into the dark sector, followed by a sinmlation of 
the dark radiation. Finally, in the last step we randomly assign the bosons to Hd, b^, ad, or 
hd and decay them accordingly, thus producing lepton jets. We take the dark fermions to be 
stable and only account for their radiation. 

3.3 Lepton jets' morphology 

Before moving on to discuss the results of the simulations, we briefly discuss the qualitative 
differences between the resulting lepton jets and the type of jets one can expect from QCD. 
Lepton jets differs from QCD jets by both composition and shape. 

The composition of the lepton jet is affected by the ratio of electrons/muons to pions 
which is determined by the decay modes of the dark vector-boson. The decay branching ratio 
of the vector-boson b^ into pions is dictated by the electromagnetic form- factor at = ml, 
also known as the R ratio [19] . If the vector-boson mass is very close to the p-meson resonance, 
then its decay is mostly into pions and speaking of "lepton jets" is not very appropriate. In 
general, however, a sizable branching fraction into leptons can be expected and in what 
follows we consider several benchmarks with different branching fractions (Br (6^ tt+tt") = 

''in what follows wc mention leptons only for the sake of brevity, but pions are also possible. We discuss 
this issue more carefully below. 
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1/7,1/3,3/5 and taking the muon vs. electron branching fraction to be equah). This is an 
important effect to model because pion contamination will reduce the efficiency for lepton jet 
searches with hadronic isolation cuts as we discuss in the next section. 

Regarding the shape of lepton jet, they are usually made of a "hard core" of tightly packed 
high energy leptons coming from the primary dark bosons, and a "soft-shell" of somewhat 
sparse and less energetic lepton pairs coming from dark radiation. This configuration is quite 
different from usual QCD jets and can be used to search for this objects as we discuss in 
the next section. The "soft-shell" itself may contain additional discriminating power. Since 
the dark photon is relatively light < GeV, even the soft radiation usually results in fairly 
collimated leptons. We therefore expect the pattern of energy deposition in the "soft-shell" 
to be of fairly isolated hits. 

It is important to realize that there are certain regions of parameter space where the 
resulting jets look very similar to QCD jets and their discovery is difficult at best. As 
we just mentioned, when the dark vector-boson mass is very close to the p-meson mass, 
it decays mostly into pions. If the dark coupling is very large, radiation is a substantial 
effect which will smear out the distinction between the "hard core" and "soft shell". In this 
case, the resulting jets resemble QCD jets in both composition (mostly mesons) and shape 
(smeared energy deposition). That said, in most other parts of parameter space, lepton jets 
are sufficiently different from QCD jets by both composition and shape that their discovery 
should be possible. 
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Figure 8: The different decay modes associated with the dark bosons: a) direct decay into leptons 
tlirough kinetic mixing; b) Hd decays into two on-shell which then decay into two leptons each; c) 
ad decays into hd and dileptons through an on-shell or off-shell dark vector-boson depending on the 
detailed spectrum; d) hd if Ughter than the other bosons typically decays outside the detector and 
constitutes missing energy. 

4. Lepton Jets and Experimental Searches 

In this penultimate section we present the results of simulations for the different processes 
discussed above. Production of the dark states was simulated using Madgraph [20]. The later 
cascade and decay back into the SM was simulated with a private code using Mathematica. 

4.1 Lepton jets from rare decays 

As discussed in the previous section, the effects of dark radiation in the case of rare decays 
into the dark sector are mild with possibly one radiated dark photon in a fraction of the events 
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(left pane of Fig. 6). The structure of the lepton jets in this case is therefore straightforward 
to understand. If the decay is into dark fermions, the event contains large amounts of missing 
energy and possibly one or two dileptons from the fermions cascade. On the other hand 
the decay into dark bosons lead to 4 distinct topologies. Assuming CP conservation, the 
must decay into one CP-even and CP-odd boson. On one side of the event we can expect 
either 6^ or and on the other side or ha- Their decays are depicted in Fig. 8 and 
the corresponding lepton jets are comprised of 2 or 4 collimated leptons with some of the 
events containing missing energy. Radiation will increase the lepton multiplicity in some of 
the events, but should not substantially affect the event topology and lepton jet structure. 
We consider this an important channel for lepton jet searches since it is rather clean with a 
simple event topology as depicted on the left of Fig. 9. 



Figure 9: The event topologies discussed in section 4. On the left we depict a rare decay into dark 
states which subsequently decay as in Fig. 8. This would result in two isolated lepton jets recoiling 
against each other. On the right is a neutralino pair production event with the neutralino ultimately 
decaying into the dark sector. We allow the dark bosons to decay as usual, but keep the dark fermions 
stable, aside from possible radiation. The events therefore consist of 2 hard lepton jets, missing energy, 
and softer leptons coming from radiation of the dark fermions (radiation from the dark bosons would 
normally be clumped together with the harder leptons coming from the cascade.). 

4.2 Lepton jets from neutralino cascades 

In the case of neutralino cascades, shown on the right of Fig. 9, the effects of dark radiation 
are more substantial. If not for radiation, one would expect two clean lepton jets in almost 
every event, coming from the scalars decaying into leptons (with the exception of which 
would constitute missing energy, and any decays involving pions). Including radiation, the 
dark fermions may also produce lepton jets (albeit softer) which would increase the lepton 
jet multiplicity. On the other hand, dark radiation might also deteriorate the signal. In 
particular, a justified concern is the pollution of lepton jets by pions coming from radiated 
dark photons. Such contaminants will make the distinction between lepton jets and QCD-jets 
difficult. 
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In order to investigate these effects in detail, we simulated the entire process, starting 
from production, going through the dark states' radiation and cascade, and ending with the 
decay back into leptons and pions. We assume that the neutralino always decays into the dark 
sector via its coupling to the dark current, Eq. (2.4). In simulating the cascades we included 
all 4 decays shown in Fig. 8 in equal proportions and kept the dark fermions stable. The results 
are presented in Tables 1 and 2 where we present the probability of finding clean lepton jets 
per event for different values of the dark gauge-coupling and the branching ratio into pions. 
Clean lepton jets are defined as follows: at least two leptons with pr > 10 GeV each in a 
cone of AR < 0.1 with hadronic and leptonic isolation of < 3 GeV in 0.1 < AR < 0.4. 

The number in brackets represents the same probability, but where leptonic isolation has 
been removed and only hadronic isolation is required. Not surprisingly, the efficiency without 
lepton isolation is higher compared with the case where isolation is required. This is especially 
so when aa is large and more radiation is expected in the isolation annulus. This effect may 
be important when looking for light sectors which are strongly coupled and radiate copiously. 
On the left of Fig. 10 we show the differential lepton jet number per event as a function of 
the lepton jet p^. As the dark coupling increases, the px tends to decrease since the stronger 
radiation dilutes the hard core of the lepton jet. On the right of Fig. 10 we depict the missing 
energy distribution in these events. It depends strongly on the mass of the neutralino and 
only weakly on the dark gauge coupling. 
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Table 1: Clean lepton jet efHciencies for different values of the dark gauge-coupling and Br (6 — )• 
tt+tt"). The neutralino mass was set to M = 100 GeV. For = dark radiation was switched off. 
The number of lepton jets increases with ad as radiation becomes more likely. The requirement for 
"clean" lepton jets, as described in the text, results in a decrease in efficiency with the growth of the 
branching ratio into pions. In brackets are efficiencies for the case where only hadronic isolation is 
required in the 0.1 < AR < 0.4 annulus. The statistical error on the efficiencies is ±0.03 

The probability of finding lepton jets decreases with increasing dark gauge-coupling since 
more radiation results in a lower average px and softer leptons. A larger branching ratio into 
pions tends to pollute the lepton jets. At Br(6 — >■ tt+tt") = 3/5, for example, it becomes very 
unlikely to observe more than one lepton jet per event. The probability of observing more 
than two lepton jets is rather small throughout the parameter space and is not depicted. 
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Table 2: Same as Table 1, but with M = 300 GeV 
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Figure 10: On the left pane we plot the differential lepton jet number per event as a function of 
their px for different values of the neutralino mass and gauge coupling. Increasing the mass tends to 
increase the average pt as expected, while a larger gauge-coupling tends to soften it because more 
radiation is generated. On the right pane we show the missing energy distribution in the events. The 
distributions are fairly insensitive to the dark gauge coupling. 



However, it is important to realize that this is to some extent a consequence of the strict 
definition of clean lepton jets. Such a strict definition is probably necessary to trigger and 
search inclusively for the harder lepton jets, but it may be desirable to relax the requirements 
somewhat for the other lepton jets in the events. In Fig. 11 we plot the probability of finding 
more than one lepton jet as a function of the pT cut on the second hardest lepton in the jet. 

Another important characteristic of a lepton jet is its lepton multiplicity. Higher mul- 
tiplicity helps with background rejection since few standard model processes can give more 
than 2 hard and collimated leptons. In Tables 3 and 4 we show the lepton multiplicity dis- 
tribution in the hardest jet for different values of the gauge-coupling and branching fraction 
into pions. A sizable fraction of the events (~ 0(10%)) have at least one lepton jet with 4 
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Number of isolated lepton jets P''"cut(^^^) 



Figure 11: On the left pane we show the dependence of the lepton jet efficiency on the isolation cut 
J2pt- On the right pane we plot the efficiency for 2 (black-solid), 3 (blue-dashed), and 4 (purple- 
dotted) clean lepton jets as a function of the PTout of the second hardest lepton in the jet. The 
parameters used are M = 300 GeV, ad = 0.1, and Br(& — >• tt+tt") — 1/7. This plot shows that by 
lowering the pt requirement on the second lepton in the lepton jets, one can increase the efficiency for 
their observation. 





Lepton 


Multiplicity 


in Clean Lepton Jets 








2 Leptons 


4 Leptons 


6 Leptons 




1/7 


1/3 


3/5 


1/7 


1/3 


3/5 


1/7 


1/3 


3/5 







0.41 


0.31 


0.21 


0.22 


0.13 


0.06 


0. 


0. 


0. 


0.01 


0.37 


0.33 


0.23 


0.23 


0.15 


0.05 


0. 


0. 


0. 


0.03 


0.33 


0.31 


0.22 


0.2 


0.14 


0.06 


0.01 


0. 


0. 


0.1 


0.29 


0.26 


0.18 


0.18 


0.12 


0.04 


0.02 


0.01 


0. 


0.3 


0.22 


0.19 


0.14 


0.09 


0.08 


0.03 


0.02 


0.01 


0. 



Table 3: Lepton multiplicity in the hardest lepton jet for different values of the dark gauge-coupling 
and Br(6 tt+tt"). The neutralino mass was set to M = 100 GcV. Odd number of leptons is fairly 
unlikely since the vector-boson always decays into a pair of leptons. 

hard leptons in it coming mostly from Hd decay. The number of such lepton jets diminishes 
as the gauge-coupling increases because dark radiation pollutes the annulus 0.1 < AR < 0.4. 
This reduction is less significant if one removes the requirement of leptonic isolation in the 
0.1 < AR < 0.4 region. In Fig. 12 we show the differential distribution of lepton jets with a 
given number of leptons as a function of pT- 
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Table 4: Same as Table 3, but for M = 300 GeV. 
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Figure 12: The average number of lepton jets with a given number of leptons in them ni^p = 1, 2, 3, 
and 4 for Mf^^ = 100 GeV (300 GeV) on the left (right) with Ud = 0.1 and Br(& ^ tt+tt") = 1/3. 
Most lepton jets contain 2 leptons. The niep = 1 case does not really satisfy the requirements needed 
to be identified as a lepton jet, but is included here to illustrate that isolated leptons in those events 
are fairly soft. They originate from soft radiated dark photons which have decayed into two, fairly 
separated leptons. This does not include the possibility of a hard isolated lepton originating from 
chargino to neutralino cascade, a possibility which was not simulated here. 



4.3 Experimental Searches 

These findings suggest that it is possible to look for lepton jets in a fairly inclusive fashion and 
compel us to put forward some suggestions for such experimental searches. One such search 
was in fact already done in the Tevatron [6] where a fairly similar definition for lepton jets as 
the one discussed below was employed. We believe that many of the conclusions arrived at 
here will remain true and carry over in the case of a non-abelian dark gauge group as long as 
it does not condense. It is useful to try and arrive at some quantitative definition of a lepton 
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jet and wc suggest the following^: 

Two or more leptons each with px > 10 GeV inside a cone of AR < 0.1 with 
hadronic/leptonic isolation cut of^pT < 3 GeV in an annulus of 0.1 < AR < 0.4 
around the lepton jet. 

We view this definition as a template to which more adequate figures can be added 
later once issues of backgrounds, triggering and other experimental difficulties are thoroughly 
analyzed and resolved. In particular, it will be important to bring the following parts under 
control: 

1. "Two or more leptons" - Muons likely suffer from less background than electrons (either 
instrumental or physical) so we may require fewer of them than we do electrons. When 
only two hard leptons are present, one might demand the hardest two leptons to have 
the same flavor and opposite sign as this may prove useful in reconstructing the tracks 
as they separate in the magnetic field®. This might be more difficult to do in the case 
of Ha decay (or non-abelian scenarios) where more than one pair of hard leptons is 
expected. 

2. "pr > 10 GeV" - This figure can be changed depending on the lepton species under 
consideration. Also, it may be refined to allow for a softer secondary lepton. 

3. "Ai? < 0.1" - This figure may be tuned depending on the process under consideration 
and the expected size of the opening angle. Also, one may want to define a lower limit 
in order to reduce background. 

4. "^Pt < 3 GeV" - A hadronic/leptonic isolation cut around the lepton jet core is 
necessary in order to fight background. It may be possible to relax the requirement for 
leptonic isolation if the background is not too large. The precise annulus in which this 
isolation cut is required (tentatively, 0.1 < AR < 0.4) should be thought of carefully 
especially since dark radiation leads to more dispersed lepton jets. 

5. "around the lepton jet" - In this work we defined the cone as follows: We first used an 
iterative algorithm to build the lepton jet 4-vector. Starting with the hardest lepton we 
collected all leptons within AR = 0.1 around it and added their 4- vectors to the lepton 
jet. This was repeated until no further leptons were found within AR = 0.1 around 
the lepton jet 4-vcctor. This same 4-vector was then used to define the isolation cone 
0.1 < AR < 0.4. This choice was motivated by the physics since the combined 4-vectors 
reconstruct the dark state. However, other definitions are possible, for example the one 
employed by Ref. [6] where several isolation cones were defined, one about each lepton 
in the lepton jet. 

®We thank the participants of Boost09 and especially, B. Demirkoz, A. Haas, R. Moore, P. Schuster, N. 
Toro, and J. Wacker for very fruitful discussions concerning these issues. 
®We thank Andy Haas for pointing this out to us. 
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6. Lepton jets can have a displaced vertex as well. This case should be carefully studied 
by itself since in general the trigger algorithms assume that electrons and muons have 
a track in the inner detector and therefore displaced lepton jets may fail to trigger^. 

This definition is likely to be revised and modified once collider effects are studied and 
a better understanding of the different backgrounds is developed. Once a proper operational 
definition of a lepton jet is formed, and maybe several interesting subclasses are identified, it 
is possible to conduct inclusive searches for these objects: 

1. Lepton jet recoiling against a QCD-jet would be an inclusive search for a prompt dark 
photon production [3]. 

2. Two lepton jets recoiling against each other and reconstructing the would be an 
interesting signal of rare decays into the dark sector and can be looked for at LEP, 
Tevatron, and LHC. 

3. Two (or more) lepton jets together with missing energy and possibly other isolated 
final states (e.g. a muon, an electron, and etc.) can be the result of electroweak-ino 
production and their eventual cascade into the dark sector. 

4. Lepton jets in association with QCD-jets could be the result of strong production of 
colored particles which eventually cascade into the dark sector. This production mode 
was not investigated in this paper, but is important to investigate in future works on 
the subject since it enjoys a very large cross-section. 

5. Conclusions 

In this paper we investigated the (supersymmetric) production of dark states through elec- 
troweak processes. We carefully included the dark radiation that accompanies highly boosted 
dark states as well as their cascade evolution back into leptons and pions. This allows for 
a more realistic characterization of the resulting lepton jets and their properties. In partic- 
ular, we put forward an operational definition for lepton jets that should aid in performing 
inclusive searches for these objects. We simulated the relevant processes (without including 
detector effects) and conclude that the efficiency for identifying clean lepton jets is rather high 
throughout the parameter space, even when possible pollution from pions and other leptons 
is included. 

The entire process depends on several disjoint assumptions and we have attempted to 
keep the discussion modular throughout. The production of dark states depends of course on 
the existence of the dark sector and its coupling to the standard model. The dark radiation 
is a result of the ^ GeV scale gauge group and depends on the gauge coupling and (logarith- 
mically) on the energy carried by the dark states. The cascade back into the standard model 

'^We thank Roger Moore for pointing this out to us. 
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depends on the details and spectrum of the dark sector which in this work we took to consist 
of two dark chiral multiplets (scalars and fermions). Finally, the decay branching fraction into 
pions vs. leptons is determined by the electromagnetic form-factor at the dark gauge-boson 
mass and we have modeled this effect by considering several pion branching fractions. 

However, while our analysis depends on numerous assumptions, wc believe that the re- 
sulting phenomenology is rather robust and depends mostly on the existence of such a light 
sector with a small, but not negligible, interaction with the standard model. As such, the 
present work should aid future searches for lepton jets in high energy colliders irrespective of 
the precise details of the light sector. 
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A. A Simple Abelian Model 

Throughout this work, we have been assuming that the WIMP is coupled to some abelian 
gauge group U{l)d which is broken and becomes massive at ~ GeV. In this appendix we 
briefly describe the dark higgs sector responsible for the spontaneous breaking of the gauge 
symmetry. Our notation follows that of Ref. [21] and the vertices can be read off the MSSM 
vertices by restriction of the electroweak higgs doublets to the neutral higgses' interactions 
with the vector-boson. Wc do not discuss issues of naturalness or other model building 
concerns [22, 23, 24], but simply use the model below as a framework to describe the resulting 
phenomenology. 

The dark higgs scalar potential, including soft terms is given by. 



V{hi,h2) = 




(A-1) 



where the soft masses include possible contributions from /x-like terms. We can express these 
in terms of the mass eigenstates, hj, and which are the light and heavy real scalars, a^, 
which is the pseudo-scalar, and bd which is the Goldstone eaten by the gauge field b^, 

hi = vi -\ — -1= {Hd cos a — hd sin a + iud sin /3 — ibd cos /3) (A-2) 
v2 

h2 = V2-\ — 1= {Hd sin a + hd cos a + iad cos P — ibd sin P) (-^"3) 
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where f 1^2 are the vacuum expectation values (VEV) that spontaneously break U (l)d, and 
their ratio determines P through tan^ = V2/V1. The angle a is determined through the 
mass matrix of the real scalars. It is convenient to adopt tan/3, 171^2, and as independent 
parameters. We then have, 

ml^ = (tan P + cot P) (A-4) 

^Ha,ha = \ (rnl^ +ml± ^J {ml^ + mlf - ^rnlml^ cos'^ 2/3^ (A-5) 

(2 2 \ 
T- (A-6) 

Most importantly, the following hierarchy is always maintained m^^ < rua^ , < , 
which strongly constrains the type of decays involved. The relation between rua^ and mi, is 
constrained by the requirement that 7715 decays into leptons and not into hd and Ud, so Ud 
cannot be taken too light compared with 7715. Therefore, in large regions of parameter space, 
the decays in the dark higgs sector are as shown in Fig. 8. 

When simulating the cascades in the dark sector we assumed equal branching ratio of 
the MSSM neutralino into the different dark states, which is a very good approximation 
considering the large hierarchy of scales. We used phase-space distributions for the decay 
and chose niH = 1-1 GeV, nib = 0.5 GeV and iria = 1 GeV for concreteness, however, the 
kinematics is determined mostly by the large energy scale involved, namely, the neutralino 
mass. 
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